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ABSTRACT 

We investigate the role of ram pressure stripping in the Virgo cluster using N-body simulations. 
Radial orbits within the Virgo cluster's gravitational potential are modeled and analyzed with 
respect to ram pressure stripping. The N-body model consists of 10 000 gas cloud complexes which 
O \ can have inelastic collisions. Ram pressure is modeled as an additional acceleration on the clouds 

4— » ' located at the surface of the gas distribution in the direction of the galaxy's motion within the 

^ \ cluster. We made several simulations changing the orbital parameters in order to recover different 

stripping scenarios using realistic temporal ram pressure profiles. We investigate systematically 
the influence of the inclination angle between the disk and the orbital plane of the galaxy on the 
^ ' gas dynamics. We show that ram pressure can lead to a temporary increase of the central gas 

H \ surface density. In some cases a considerable part of the total atomic gas mass (several 10* Mq) 

can fall back onto the galactic disk after the stripping event. A quantitative relation between the 
orbit parameters and the resulting Hi deficiency is derived containing explicitly the inclination 
angle between the disk and the orbital plane. The comparison between existing Hi observations 
and the results of our simulations shows that the Hi deficiency depends strongly on galaxy orbits. 
It is concluded that the scenario where ram pressure stripping is responsible for the observed Hi 
deficiency is consistent with all Hi 21cm observations in the Virgo cluster. 

Subject headings: ISM: clouds - ISM: kinematics and dynamics - Galaxies: cluster: individual: Virgo 
Cluster - Galaxies: evolution - Galaxies: interactions - Galaxies: ISM - Galaxies: kinematics and 
dynamics 

1. Introduction 



^also at: Max-Planck-Institut fiir Radioastronomie, Auf 

dem Hugel 69, D-53121 Bonn, Germany Since Gunn & Gott (1972) have introduced the 

concept of ram pressure stripping, which can affect 
galaxies moving inside the Intracluster Medium 
(ICM) of a galaxy cluster, this mechanism has 



been invoked to explain different observational 
phenomena as the HI deficiency of spiral galaxies 
in clusters (Chamaraux et al. 1980 , Bothun et al. 
1982, Giovanelli & Haynes 1985) or the lower star 
formation activity of cluster spiral galaxies (e.g. 
Dressier et al. 1999, Poggianti et al. 1999). 

Some other mechanisms have been proposed af- 
fecting spiral galaxies in clusters, among which: 

• galaxy-galaxy interactions (galaxy "harass- 
ment"; Moore et al. 1996, 1998) 

• galaxy-cluster gravitational interaction (Byrd 
& Vahonen 1990, Valluri 1993) 

All these kinds of interaction have in common that 
they can in principle remove the neutral gas from 
a galaxy inducing a decrease of its star formation 
activity. 

The strength of these interactions depends cru- 
cially on galaxy orbits. Radial orbits allow galax- 
ies to go deeper into the cluster potential where 
their velocity increases considerably and where the 
galaxy density and the density of the ICM is sub- 
stantially higher. Dressier (1986) analyzed Gio- 
vanelli & Haynes (1985) data showing that Hi de- 
ficient spiral galaxies are on radial orbits. Numer- 
ical simulations (Ghigna et al. 1998) showed that 
galaxy halos evolving within the cluster settle on 
isotropic orbits with a median ratio of pericentric 
to apocentric radii of 1:6. 

The recent analysis of a spectroscopic catalog 
of galaxies in ten distant clusters (Dressier et al. 
1999, Poggianti et al. 1999) has shown that the 
galaxy populations of these clusters are charac- 
terized by the presence of a large number of post 
starburst galaxies. Poggianti et al. (1999) con- 
cluded that the most evident effect due to the 
cluster environment is the quenching of star for- 
mation rather than its enhancement. They found 
two different galaxy evolution timescales in clus- 
ters, (i) A rapid removal of the gas (^ 1 Gyr). 
(ii) A slow transformation of morphology (several 
Gyr). They stated that the mechanism responsi- 
ble for the fast gas removal without changing the 
galaxy's morphology is very probably ram pressure 
stripping. Furthermore, Couch et al. (1998), using 
HST data, claimed that there are fewer disturbed 
galaxies than predicted by the galaxy harassment 
model. They suggested that ram pressure strip- 
ping by the hot ICM truncates star formation. 



In the local universe, the existence of a large 
number of post starburst galaxies in the Coma 
cluster seems to indicate a triggering of the star 
formation activity by effects which are related to 
the cluster environment (Bothun & Dressier 1986, 
Caldwell et al. 1993). Caldwell & Rose (1997) 
observed an enhanced number of post starburst 
galaxies within cluster environments. On the con- 
trary, Zabludoff et al. (1996) found that interac- 
tions with the cluster environment, in the form of 
the intracluster medium or cluster potential, are 
not essential for "E-l-A" formation. HST ivaages 
of the Coma starburst and post starburst galax- 
ies (Caldwell, Rose, & Dendy 1999) revealed that 
star formation takes mainly place in the inner disk 
of these galaxies and is suppressed in the outer 
disk. If the starburst is related to the cluster en- 
vironment, ram pressure stripping can naturally 
account for these characteristics. 

The best place to study the gas removal due 
to ram pressure is the Virgo cluster as it is the 
closest cluster which can be observed in great de- 
tail. Most of the spiral galaxies seem to have 
entered the cluster only recently (within several 
Gyr, TuUy & Shaya 1984). About half of them 
became Hi deficient (Giovanelli & Haynes 1985). 
Their Hi disk sizes are considerably reduced (van 
Gorkom & Kotanyi 1985, Warmels 1988, Cayatte 
et al. 1990, 1994). Concerning the stellar con- 
tent, their intrinsic color indices arc not signif- 
icantly different from field galaxies of the same 
morphological type (Gavazzi et al. 1991, Gavazzi 
et al. 1998). But despite the Hi deficiency, cluster 
galaxies do not show a reduced CO content (Ken- 
ney & Young 1989, Boselh et al. 1997) neither 
a reduced infrared luminosity (Bicay & Giovanelli 
1987). This means that the neutral gas was re- 
moved without affecting the molecular component 
and without changing colors very much. The only 
interaction cited above producing this signature is 
ram pressure stripping. 

Very few simulations have been done to quan- 
tify ram pressure stripping using Eulcrian hydro- 
dynamic (Takeda, Nulsen, & Fabian 1984, Gaetz, 
Salpeter & Shaviv 1987; Balsara, Livio, & O'Dea 
1994) or SPH codes (Tosa 1994; Abadi, Moore, 
& Bower 1999). We use a sticky particles model 
in which each particle represents a cloud complex 
with a given mass dependent radius. The viscos- 
ity of the clumpy ISM is due to inelastic collisions 



between the particles. The effect of ram pressure 
is modeled as an additional acceleration on the 
particles located in the direction of the galaxy's 
motion. 

This article is devoted to the question if ram 
pressure stripping can account for the observed 
Hi deficiency, Hi distributions, and velocity fields 
of spiral galaxies in the Virgo cluster. Further- 
more we investigate possible star formation mech- 
anisms related to a ram pressure stripping event. 
We first discuss radial galactic orbits within the 
cluster gravitational potential (Section 2). The 
model simulating the neutral gas content of an in- 
falling galaxy is described in Section 3, followed 
by the detailed description of the simulations (Sec- 
tion 4). The evolution of the stripped gas in the 
intracluster medium is investigated in Section 5. 
In Section 6 we give the results of these simula- 
tions and discuss them in Section 7. We confront 
observations of various quantities of spiral galax- 
ies with the knowledge acquired with the help of 
the simulations (Section 8) and give possible ex- 
planations for their dynamical and physical state 
(Section 9 and 10). The summary and conclusions 
are given in Section 11. 

We adopt a velocity of 1150 kms""'^ for the 
Virgo cluster (Huchra 1988) and a distance of 
17 Mpc. 

2. Galaxy orbits 

The hot intracluster medium (ICM) of density 
piCM exerts a ram pressure: Pram = Picm ^'laiaxy 
on a galaxy, which moves with a velocity iigaiaxy 
through it. With decreasing distance to the clus- 
ter center the galaxy's velocity, the ICM density, 
and thus ram pressure increase. Therefore, the ef- 
ficiency of ram pressure stripping depends strongly 
on the shape of the galaxy orbit within the cluster. 
Radial orbits lead a galaxy deeper into the clus- 
ter core where ram pressure stripping is the most 
important. We investigate in this section the influ- 
ence of the shape of radial orbits on the efficiency 
of ram pressure stripping. 

The gas in the galactic disk is removed if the 
ram pressure of the ISM is greater than the restor- 
ing gravitational force per unit area provided by 
the galaxy's disk (Gunn & Gott 1972). The strip- 
ping radius _Rstr can be estimated using the equa- 
tion established by Gunn & Gott (1972) with S, — 
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R-^ (Binney & Tremaine 1987): 



Sgas Vrot^ — PICMW. 



galaxy : 



(1) 



where Yj, 



/gas 



is the stellar/gas surface den- 



sity, picM is the ICM density, -yrot is the ro- 
tation velocity, and Wgaiaxy is the galaxy's ve- 
locity within the cluster. We assume Sgas = 
const = 10^^ cm^^ and v^ot = 150 kms~^. 
The minimum ram pressure p™^ to affect the 
atomic gas disk can be determined by setting 
R = i?Hi = 15 kpc. For our model galaxy this 
leads to vfl^ ^ 500 cm'^ (kms-i)^. Fig. 1 
shows the galaxy velocity as a function of the 
ICM density for constant values of the ram pres- 
sure Pram = PiCM^gaiaxy ^hc area bclow the 
solid line corresponds to the part of the parame- 
ter space where PiCMWgaiaxy <500 cm^^ (kms~^)^, 
i.e. where ram pressure is not effective. The curves 
correspond to a values of picM^gaiaxy— 1000' 2000, 
5000, 10000 cm^'^ (kms^^)^ that we used in our 
simulations. For the simulations we varied the in- 
clination angle O between the disk and the orbital 
plane: 6=0°, 20°, 45°, 90°. With this defini- 
tion of the inclination angle, 0=0° means edge-on 
stripping and 6=90° face-on stripping. For our 
model, this represents the physically most mean- 
ingful definition and has the advantage that it 
leads to simple fitting functions of our numerical 
results. 

We took a /3-model (Cavaliere & Fusco-Femiano 
1976) to describe the total mass density and the 
ICM density: 



/'oll + ir) 



(2) 



where tq is the core radius, /3 is the slope pa- 
rameter, and po is the central density. The val- 
ues for the ICM are: /3=0.5, rc = 13.4 kpc, and 
po = 4 X 10~^ cm^'^. For the total mass profile we 
used /3=1, rc=0.32 Mpc, and pc = 3.76 x 10-* 
Mqpc^^. This gives a total mass of Mtot = 
1.4 X 10^^^ M© at a radius of 1 Mpc (Schindlcr, 
Bingelli, & Bohringer 1999). We have calculated 
orbits of test particles within the gravitational po- 
tential given by this mass distribution. 

The giant elliptical galaxy M86 has a total 
mass of about one tenth of the whole cluster mass 
(Schindler et al. 1999). The perturbation of the 
total gravitational field due to M86 is thus not neg- 
ligible. We have added the gravitational potential 



corresponding to M86 (/3=0.5, re— 74.2 kpc, and 
PC = 6.26 X 10"'^ Mopc"'^). First, we have made 
test runs where we have put M86 on an eUipti- 
cal orbit around M87. The resuhs for the orbit- 
ing perturber (M86) are that (i) M87 begins also 
to move and (ii) the additional mass distribution 
perturbs the trajectories of the test galaxy pro- 
ducing a wide variety of radial orbits with differ- 
ent impact parameters. Without the perturbing 
gravitational potential the orbits would lie in one 
single plane and the impact parameter, i.e. the 
closest approach to the cluster center, would stay 
constant. Second, we have fixed M86 in space. In 
this case the galaxy orbits are perturbed in the 
same way as in the previous case. We can thus 
conclude that in both cases the presence of M86 
near M87 changes significantly the trajectories of 
galaxies on radial orbits. For simplicity we have 
chosen to fix M86 in space. 

Fig. 2 shows the orbits calculated with the 
fixed gravitational potential of M87 and M86. 
The asymmetric gravitational potential leads to 
a variation of the impact parameter and the max- 
imum velocity for each orbit. The initial con- 
ditions were fixed such that the resulting max- 
imum ram pressure varies between 1000 and 
5000 cm~'^ (kms~^)^. If a spiral galaxy is grav- 
itationally bound to the cluster, these high ec- 
centricities are needed in order to obtain the 
highest ram pressure maxima. However, for val- 
ues of the maximum ram pressure smaller than 
~3000 cm^'^ (kms~^)^ slightly less eccentric orbits 
are also possible. The smallest impact parameter, 
i.e. the closest approach to the cluster center, is 
~90 kpc with a maximum velocity of 2200 kms^^. 
A small change of the initial conditions, resulting 
in an impact parameter of ~70 kpc, leads to a 
maximum ram pressure of 10000 cm~^ (kms^^)^. 
Given the uncertainties of the determination of 
the enclosed mass and the ICM density profile, we 
think that this extreme orbit is still realistic for 
the Virgo cluster. The largest impact parameter 
is ~200 kpc. The timescale between two passages 
through the cluster core is about 10^° yr. A galaxy 
will therefore have at most a few passages during 
its lifetime within the cluster. Thus the simulation 
does not trace the trajectory of one single galaxy 
but creates a variety of orbits corresponding to 
different initial conditions. Nevertheless, it is not 
excluded that a galaxy can have two passages 



through the cluster core. 

We define the normalized ram pressure as 

Pnorm = (PlCM^gaiaxy )/(PoWo ) , (3) 

where po=10"^ cm^^ and uo=1000kms^^ are rep- 
resentative averaged values. The resulting normal- 
ized ram pressure of the orbits shown in Fig. 2 can 
be seen in Fig. 3. This graph illustrates that the 
perturbation of the galaxy orbits by M86 can re- 
sult in a change of the ram pressure strength up 
to 50% because of the small core radius of the 
ICM distribution. A small change of the impact 
parameter leads to a considerable change of the 
ICM density, whereas for values of the impact pa- 
rameters between 100 and 200 kpc the change of 
the maximum galaxy velocities is only of the or- 
der of ^10%. For these radial orbits ram pres- 
sure dominates the thermal pressure of the ICM 
(T ^ 10^ K) for distances less than one Abell ra- 
dius (~ 2.5 Mpc) as can be seen in Fig. 4. 

In order to compare these orbits directly with 
observations, we plot in Fig. 5 the line-of-sight 
(LOS) velocity with respect to the cluster mean 
velocity as a function of the projected distance to 
the cluster center (M87). The number of orbits is 
sufficient to be representative for radial orbits over 
the whole volume of the Virgo cluster. The region 
above the envelope of these curves corresponds to 
more circular orbits. This diagram will be used 
later on to compare the model orbits directly with 
observations. 

We can now address the question if we can 
constrain the orbital shape of Hi deficient galax- 
ies. As stated above, a minimum ram pressure of 

P^m = PiCMWgalaxy ~ 500 Cm"'"' (kmS-l)^ is nCC- 

essary in order to affect the atomic gas disk. In 
the case of circular orbits the galaxy velocity has 
the Keplerian value Ugaiaxy = ^k = ^i?(d$/di?), 
where $ is the gravitational potential of the clus- 
ter. This means that in the Virgo cluster a galaxy 
with a rotation velocity of Wrot ^ 150 kms~^, 
which is on a circular orbit, must be closer than 
250 kpc (50') to the cluster center in order to be 
affected significantly by ram pressure. If spiral 
galaxies have been accreted only recently, it is not 
probable that they are on circular orbits with a 
such small distance to the cluster center. On the 
observational side, Hoffman, Hclou, & Salpeter 
(1988) found a few Hi rich dwarf galaxies within 
a radius of 5° around M87. These galaxies are 



much more vulnerable to ram pressure stripping, 
because of their lower gravitational potential. At 
least some of these dwarfs must evolve on circular 
orbits with radii smaller than 1 Mpc, where they 
have never been stripped by ram pressure. 

In the case of highly eccentric orbits, the maxi- 
mum orbital velocity is given by the escape veloc- 



ity Wgalaxy ~ V2wk 



The impact param- 



eter must be smaller than 400 kpc (1.3°) in order 
to obtain Pram > 500 cm^'' (kms^^)^. For more 
massive galaxies (wrot > 150 kms~^) the maxi- 
mum impact parameter is still smaller. The mini- 
mum impact parameter of 90 kpc excludes a grav- 
itational interaction of the galaxy with M87. 

The impact parameters of the orbits shown in 
Fig. 2 are in the range between 90 and 200 kpc 
giving rise to maximum ram pressure Pram =1000- 
5000 cm~'^ (kms^^)^. We will show that galaxies 
on these orbits will have an Hi deficiency DEF > 
0.3 after a passage through the cluster center. 

3. The model 

There are different ways to model the interstel- 
lar gas numerically: (i) hydrodynamical, (ii) SPH, 
and (iii) sticky particle models. Each class has its 
advantages and limits. Since the ISM is neither 
a continuous medium (hydrodynamics, SPH) nor 
exclusively made of clouds (sticky particles), one 
has to chose the class which is well adapted for the 
investigated astrophysical problem. In the follow- 
ing we will motivate the choice of our model. 

The interstellar Hi gas is a clumpy medium (see 
e.g. Kim et al. 1998 for the LMC). In the Galaxy 
roughly half of the interstellar Hi is in the warm 
phase (r ^ 6000 K) with a volume filling factor 
of ^0.5. The other half is in the form of cold 
(T ~80 K) and dense clouds with a filling fac- 
tor of ~0.01 (Kulkarni & Hciles 1988). If they 
are in pressure equilibrium with the rest of the 
ISM at nT = 3000 Kcm"^, they have a density 
of n ~40 cm^''. The largest clouds have column 
densities of ^ 2 x 10^° cm^^. On the largest scale, 
three dimensional structures of the atomic gas 
phase are major Hi concentrations around dark 
clouds and newly formed stars. They have col- 
umn densities of iVni ^ 10^^ cm^'^ and masses of 
~ 10^ Mq (Kulkarni & Heiles 1988). Thus, the 
cold atomic phase of the ISM is well represented 
by discrete particles. Concerning the warm neu- 



tral phase, Murakami & Babul (1999) have shown 
in their hydrodynamical simulations of the effects 
of ram pressure on an elliptical galaxy that the 
continuous ISM tends to be stripped in the form 
of distinct clouds if they included the external ICM 
pressure. These clouds have densities n ~ 1 cni^"^ 
and radii Rd ^ 300 pc, resulting in a column den- 
sity of A^Hi ~ 10^^ cm""^. We therefore think that 
the warm neutral phase, which is subjected to an 
external pressure due to the surrounding ICM and 
to ram pressure forces, can also be modeled by dis- 
tinct particles. 

These considerations motivate our approach to 
model the entire atomic phase of the ISM as sticky 
particles where each particle represents a gas cloud 
complex. We assume that the clouds have a con- 
stant column density. The relevant physical prop- 
erty of the ISM with respect to ram pressure is the 
column density of the gas/clouds. It is adapted to 
observed values. 

The effect of ram pressure is simulated as an ad- 
ditional acceleration applied on the clouds which 
are exposed to the ICM in the direction of the 
galaxy's motion within the cluster. 

3.1. The model galaxy 

Since we do not want to include gravitational 
interactions, we let the gas clouds evolve in a given 
analytical gravitational potential. It consists of 
two spherical components and a disk (Allen & San- 
tillan 1991): 

• the dark matter halo potential which is given 

by 

R 1.02ahaio 



1.02 



L l + (i?/ahalo) 

where 



^^+ln(l+(i?/ahaio)'-°')] 



l.02^lloo 



M(R) = - 



^halo(-fi'/ahalo) 



+ (i?/ahaio)i-02 



(5) 



(6) 



R is the distance to the galaxy center, 
Mhaio = 8.6 X lOi" Mq, and ahaio=12 kpc. 



• The bulge potential which is given by 

^bulgo 



^bulgo(-R) — — 



-^ + "bulge 



(7) 



where Mbuigo — 5.6x10'' M0 and 6buigo=387 
pc. 

• The disk potential which is given by 



*disk(?', z) 



^r^ + (fldisk + ^z^ + 6^1,5" ' 

(8) 
where R^ = r^ + z"^ , Afdisk=2.6xl0i" Mq, 
adisk=2.7 kpc, and fedisk=250 pc. 

We fix the optical diameter of the disk at i'opt =20 kpc 
(^7.5 disk scale lengths). We have included self- 
gravity between the cloud complexes. Fig. 6 shows 
the resulting rotation curves after a disk evolution 
of 10^ yr for the different components as well as 
the decomposition and the total model rotation 
curve. The latter is used as initial conditions for 
our simulations. We thus obtained a constant 
rotation curve of ^140 kms^^. This rotation ve- 
locity corresponds to a medium mass spiral galaxy. 
Fig. 7 shows the Hi cloud distribution in z direc- 
tion after 10^ yr. We end up with a realistic height 
of the Hi layer of ^^1000 pc at a radius of 10 kpc. 
The initial Hi diameter is Dj^j"=30 kpc. 

The cloud mass distribution is n{m) ex tti^^'^ 
(Knude 1981). The initial total gas mass of the 
galaxy is 6 x 10^ M0, the total number of clouds 
is 10 000. We adopt a constant Hi column density 
for all clouds (Sanders, Scoville & Solomon 1985) 
which is fixed at Sci = 7.5 x 10^° cm~^. This value 
is consistent with the measurements made by Rots 
et al. (1990) on nearby galaxies. Assuming the 
surface density to be constant has the advantage 
that the acceleration due to ram pressure is the 
same for all clouds disregarding their masses. This 
leads to the following mass-radius relation which 
is relevant for the collision rate: 



^coU 
cloud 



V^cioud/(7rAE) 



(9) 



where E is the gas surface density of one cloud. 
In reality the number of Hi clouds is much larger 
than the number of particles in our model. To 
compensate this effect we have chosen the factor 
A in the way that the mean free path of one cloud 
is the same for 10 000 clouds as for a realistic cloud 
number. 

In addition, in the interior of the galaxy gas 
clouds become denser building molecular cores. 
This results in smaller cloud sizes and higher sur- 
face densities. As the atomic gas represents the 



outer layer of these clouds, they can be hardly 
pushed by the ram pressure. We take this effect 
into account by adding a factor 7 in the equation 
for the acceleration due to ram pressure: 

aram = Pram/il Ed) , 7=15 CXp ( - (r/ro)) + 1 

(10) 
where ro—2 kpc (~ 0.4" for the Virgo cluster). 
This profile is consistent with the CO (1-0) obser- 
vations of Virgo cluster spirals (Kcnney & Young 
1988). 

When orbiting around the galaxy center, the 
clouds can have inelastic collisions. During these 
collisions they can exchange mass and larger 
clouds can grow through coalescence. The out- 
come of a collision can be one (coalescence), two 
(mass exchange), or three fragments (fragmenta- 
tion). Let the radius of the first cloud be ri, that 
of the second cloud r2- Let the impact parameter 
be b, the velocity of the fragment Vf , and the es- 
cape velocity Ucsc- We follow the prescriptions of 
Wiegel (1994): 

• for ri — r2 < b < ri + r2'. 
fragmentation 

• for b < ri — r2 and Wcsc > vf- 
mass exchange 

• for b < ri — r2 and Uosc < vt'- 
coalescence 



This description is similar to that of Combes & 
Gcrin (1985) or Klaric (1995). The search for 
the next neighbors is done with the help of a 
treecode (Barnes & Hut 1986). The integration of 
the ordinary differential equation is done with the 
Burlisch-Stoer method (Stocr & Burlisch 1980) 
using a Richardson extrapolation and Stoermer's 
rule. This method advances a vector of dependent 
variables y{x) from a point x to a point x4--ff by a 
sequence of n substeps. Thus, the initial timestep 
H is divided subsequently into n=2, 3, 4, etc. sub- 
steps. At the end the solution of y{x + H) is ex- 
trapolated and an error can be estimated. The 
size of the timestep is adaptive and linked to the 
estimated error of the extrapolation. This error is 
normalized by the values of the distance covered 
during the last timestep and the velocity of each 
particle. The error level for acceptance of the ex- 
trapolated solution is a free parameter and has to 
be adapted to the physical problem treated. For 



a relative error level e < 0.1 the Courant criterion 
is fulfilled for each timestep. 

The collisions are evaluated at each timestep h — 
H/n and only those which appear for all sequences 
n are taken into account. We adopt the strategy 
to chose the error level in a way to have timesteps 
of the order of several 10^ yr. 

Ram pressure is included as an additional ac- 
celeration of those clouds which are exposed to 
the ICM in the direction of the galaxy motion. 
Fig. 8 illustrates this effect: clouds located inside 
the disk are protected by clouds on the surface due 
to their finite size (this depends of course strongly 
on the inclination angle of the galaxy with respect 
to the orbital plane). As it can be seen in Fig. 2 
the galaxy trajectories near the cluster center are 
quasi linear. This means that we do not have to 
change the wind direction during the simulation. 

4. The simulations 

Fig. 9 shows the profile of two ram pressure 
stripping events of Fig. 3 as a dashed line. In 
our simulations we approximated this ram pres- 
sure profile with a Lorentzian profile (Fig. 9 solid 
line): 

Pram ^ Pmax X (t + ^Hw) ' (-'--'-) 



where Pmax is the maximum ram pressure at the 
smallest distance to the cluster center and inw is 
the duration of the event. The simulations were 
started at t=-500 Myr. For a fixed Pmax, taw is 
chosen in order to obtain a constant ram pressure 
at i=-500 Myr: pram(~500Myr) = 0.5po^'o- This 
procedure reproduces best the ram pressure pro- 
files of the model orbits (Fig. 3). 

These values together with the wind direction 
and the resulting final parameters for each run can 
be seen in Table 1. The columns are the follow- 
ing: Col. (1): RUN. Col. (2): Maximum ram 
pressure in units of po^^o • Col. (3): Inclination an- 
gle & with respect to the orbital plane in degrees. 
We want to recall here that 0=0° means edge- 
on stripping and 0=90° face-on stripping. Col. 
(4): Time between the beginning of the stripping 
(defined by the moment when the total gas mass 
within 25 kpc Mtot < 5 x 10^ Mq) and the galaxy's 
closest passage to the cluster center. Col. (5): Ra- 
dial distance corresponding to Col. (4) assuming a 
constant velocity of the galaxy with respect to the 
cluster center of fgai=1700 kms~^. Col. (6): Final 



Hi diameter in kpc. Col. (7): Final central surface 
density in M©pc~^. Col(8): Final Hi deficiency 
(see Sect. 6.1). Col. (9): Re-accreted gas mass in 
10^ M0. Col. (10): Re-accreted gas mass divided 
by the stripped gas mass. The final Hi diameter is 
defined by a limiting column density of 10^° cm^^. 
Following Cayatte et al. (1994), the radial profile 
was determined by averaging the model distribu- 
tion in rings of integrated column density. The 
central surface density is defined as the mean den- 
sity over the inner part of the disk out to a quar- 
ter of the optical radius. For the calculation of 
the stripped and the re-accreted gas mass we used 
the enclosed mass within a radius of 20 kpc and a 
height of 1 kpc. For RUN A(10,0) 2.4 x 10*^ Mq of 
atomic gas are pushed to radii greater than 20 kpc 
and then fall back within 400 Myr. Thus, the ra- 
tio between stripped and re-accreted gas is greater 
than one. 

We have made 16 runs using different disk in- 
clinations with respect to the orbital plane and 
different values for the maximum ram pressure. 
We only use values for the maximum ram pres- 
sure which are high enough to show a visible de- 
formation of the Hi distribution. In Fig. 10 we 
show the final Hi radius i?Hi = -Dhi/2 of our sim- 
ulations as a function of the linearized estimated 
radius using the formula of Gunn & Gott (1972): 

RgG = SgasWrot/(PlCMWgalaxy) + 2 kpC, whcre 

Sgas = 7.5 X 10^" cm~2 and Wrot = 140 kms^^ 
The additional length is due to the length scale 
of the molecular gas distribution (7). The solid 
line corresponds to Rqg — Rm- As expected, the 
estimated stripping radius using Gunn & Gott's 
formula is well reproduced for an inclination angle 
O = 90° (face-on stripping) . 

In order to generalize Gunn & Gott's formula 
for inclination angles Q ^ 90°, we searched for 
a simple correlation between the stripping radius 
and (sin 9, pv^). Since we did not succeed, we con- 
clude that it can be dangerous to derive the strip- 
ping parameters only on the basis of the final Hi 
diameter. In particular, taking the velocity com- 
ponent perpendicular to the disk plane Pram = P'v\ 
does not reproduce the final Hi radii of our sim- 
ulations. We will show in Section 7.4 that the Hi 
deficiency is much better suited. 

In order to give an idea of how ram pressure acts 
in detail. Fig. 11 shows an example for an edge-on 
stripping process (RUN C(50,0)). The ram pres- 



sure maximum occurs at t=0 yr. The wind di- 
rection is indicated by the arrows whose length is 
proportional to ram pressure. At t ^ —100 Myr, 
a density enhancement begins to grow in the di- 
rections of the galaxy's motion. At t ~70 Myr, 
the gas whose rotation velocity is parallel to the 
wind direction is accelerated and driven out of 
the galaxy leading to an arm-like structure par- 
allel to the wind direction. At the opposite side 
where the gas clouds are decelerated by the wind, 
a second much less pronounced arm forms. At 
i=150 Myr ram pressure has completely ceased 
and the dynamics are determined by rotation and 
re-accretion. The main arm moves to the north- 
west of the galaxy due to its initial angular mo- 
mentum. At the same time, the gas begins to fall 
back to the galaxy in the north and south-west. 
This infalling gas is forming a second arm in the 
south-west. At t ^300 Myr a density enhancement 
in the east is observed which is due to the infalling 
material from the north. It is interesting to notice 
that the gas within the disk of the galaxy forms a 
leading m=l spiral structure. 

We show snapshots of RUN J(20,45) in Fig. 12 
as an example of an out of plane stripping be- 
ginning at with time differences of At ^ 80 Myr. 
The galaxy is seen edge-on and is moving in the 
south-east direction, i.e. the wind is coming from 
the south-east. This wind direction is indicated 
by the arrows whose length is proportional to ram 
pressure. The ram pressure maximum occurs at 
i=0 yr. At t ^-90 Myr the disk is already slightly 
deformed in z-direction. When ram pressure is 
at maximum the outer parts of the gas disk are 
pushed in the wind direction. Ram pressure ex- 
ceeds the threshold of Pram = ^PqVq early (t ~- 
300 Myr) for the given Lorcntzian profile. An im- 
portant distortion of the gas distribution can be 
observed at t ~0 yr. The gas is driven out of 
the galaxy in the direction of the wind and ro- 
tates at the same time. This creates the west- 
ern arm which develops between t=-10 Myr and 
i=100 Myr. Since the duration of the ram pres- 
sure event (tnw) is a considerable fraction of the 
rotation period, the out of plane gas in the region 
where its motion due to rotation and the motion 
of the ICM are opposite is compressed. This com- 
pressed region corresponds to the arm which de- 
velops in the north of the galaxy. For t >200 Myr 
the wind has completely ceased and the evolu- 



tion of the gas distribution is due to rotation and 
re-accretion of the stripped gas. At the end of 
the simulation (f=500 Myr) a considerable part 
of the stripped ISM is located within a radius of 
15 kpc outside the disk. Even at t ~500 Myr after 
the galaxy's closest passage to the cluster center 
the column density of the gas located outside the 
galactic disk is of the order of 10^° cm~^. 

5. The fate of the stripped gas 

We now address the question if such a large 
structure can be observed and in which phase the 
stripped gas is re-accreted. Murakami & Babul 
(1999) found that, if the thermal gas pressure of 
the ICM is taken into account, massive (10^-10^ 
Mq) neutral gas clouds are stripped. They specu- 
late that they might be able to survive in the hot 
ICM. Can an Hi blob survive long enough (i.e. 
during several 10* yr) in a hostile environment 
like the ICM? In order to estimate the evapora- 
tion time of a stripped gas cloud, we follow Cowie 
& McKce (1977). The classical approximation for 
the evaporation rate breaks down when the elec- 
tron mean free path becomes comparable or larger 
than the temperature scale height. In this case the 
classical heat flux exceeds the saturated heat flux 
and 
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where Ticm is the temperature of the ICM, nicM is 
the ICM density, and Rpc is the cloud radius. As- 
suming _Rpc=10 pc, TicM=10^ K, and riicM^lO^'* 
cm"'^ leads to cro=422. The saturated evaporation 
time is given by 
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Cap = 2.8 X lO'{nJnicM)Rpc/{Tl^M 2.73a'o' ) Y^ 

(13) 
For a given cloud column density Nc this gives 
tlllp ~ 10^(7Vc/1020 cm-2) yr. With an initial col- 
umn density of A'c = 10^^ cm^^, a stripped cloud 
evaporates within 10* yr. In the case of a mag- 
netic field configuration which inhibits heat fiux 
(e.g. a tangled magnetic field), this evaporation 
time can increase significantly (Cowie, McKee, & 
Ostriker 1981) and might attain several 10* yr. 

On the other hand, if a stripped cloud leaves the 
galaxy, its main heating source, the stellar FUV 
radiation field, disappears. The thermal equilib- 
rium gas temperature is then determined by heat- 



ing due to the soft X-ray emission of the ICM and 
coohng due to the infrared hne emission of Cii 
and Oi. Thus, when the cloud is stripped, cooling 
dominates heating, the gas cools down, becomes 
denser and begins to form H2 within a tiniescale 
of TH^ ~ lOVnc yr (HoUcnbach & Tielens 1997). 
The only atomic hydrogen resides in the outer part 
of the clouds where the X-ray flux dissociates the 
molecular hydrogen. We will now estimate the 
column density of this Hi layer. 

Maloney, HoUenbach, & Tielens (1996) have 
modeled the physical and chemical state of dense 
neutral gas exposed to an external X-ray source. 
We will extrapolate their results, which they have 
obtained for higher gas densities and X-ray fluxes. 
These authors assumed a power low X-ray spec- 
trum, whereas in our calculations the radiation 
field is made of thermal bremsstrahlung. Despite 
this difference we can use their results, because for 
hot gas temperatures of order 10^ K, the normal- 
ization constant for a given X-ray flux is not very 
different from that for a spectral index a ~ 0.7 
adopted by Maloney et al. (1996). In order to ap- 
ply their calculations we have to make sure that 
photons with energies less than 100 eV do not con- 
tribute to the ionization of the molecular cloud. 
This translates to the condition that the ionized 
column density due to photons of i? < 100 eV 
must be much smaller than the total column den- 
sity. 

The ionizing photon flux is given by the integral 
of dE/E over the range from 13.6 eV to 100 eV: 

^ ln(0.1/0.0136)Fx_,,^^,,8^ „^ .. 
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= 8.3xlO*Fx photons cm 

(14) 

where Fx = (2 — 0.5)Fo is the X-ray flux in the 
ROSAT band. In equilibrium this gives rise to an 
ionized column density 
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where a = 2.6 x W-^^iT/lO"^ K)-"-^ cm^s^^ is 
the recombination coefficient and He is the electron 
density. Using Fx = 1.5 x 10"''' ergcm^^s"^ (see 
below), T = 2 X 10'^ K, and ric = 100 cm"^ gives 
A'e ^ 10^^ cm"^. This is much smaller than the 
total cloud density of ^ 10^" cm"^. 

Maloney et al. (1996) showed that the transi- 
tion between atomic and molecular hydrogen oc- 
curs if Hx/ric > 10"^® crgscm'^s"^, where Hx is 



the energy deposition rate per particle and ric is 
the local cloud density. This deposition rate is 
approximately 



Hx -- 5.8 X IQ-^^ Fx{ 
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10^^ cm 2 



ergss 



(16) 



where A^h is the hydrogen column density. In or- 
der to estimate the X-ray flux, we take the X- 
ray surface density profile given by Schindler et 
al. (1999) 



F = Fo(l + (r/reorc)')' 



(17) 



with a core radius rcoro=13.35 kpc. Fq is deter- 
mined by 27r/p°"' rS dr = Lx, where rout =1-5 Mpc 
is the outer radius and Lx = 8.3 x 10^"^ ergss"^ is 
the total X-ray luminosity (Bohringer et al. 1994). 
This gives Fq — 1.15 x 10~^ ergscm"^ s~-^. At a 
distance of r=0.15 Mpc from the cluster center the 
X-ray flux is F ~ 1.5 x 10^^ ergscm^^s"^. With a 
hydrogen density of n = 100 cm"*^ the II/II2 tran- 
sition thus occurs at a hydrogen column density of 
Ah ^ 2 X 10^" cm"^. If the neutral gas is in pres- 
sure equilibrium with the ICM (Tjcm = 10^ K, 
nicM ^ 10^'^ cm^'^ at r=0.15 Mpc), it has a tem- 
perature of T - 100 K. 

The critical mass of an isothermal sphere (see 
e.g. Spitzer 1978) with an external pressure at the 
boundary of Pcxt/fce = "^icmTicm = 10^ cm^^K 
is Merit — 1-5 X 10^ M0. If the magnetic field 
is of the order of a i? ~ 5/iG, the critical mass 
might be as large as Merit ^ 5 x 10^ Mq. If the 
temperature of the neutral gas is r=200 K this 
critical mass becomes Merit ^ 10"* Mq. 

Such a cloud has a size of ^ 10 pc and a total 
column density of A'tot ^ 3 x 10^^ cm"^. Thus, 
only 10% of the total column density of these 
clouds is in atomic form visible in the Hi 21cm 
line. The total column density of the stripped 
gas in Fig. 12 observed with a 20" beam (VLA 
C-l-D array configuration at 21 cm) is Aj^'^"'^''" ^ 
several 10^° cm"^. Without evaporation the ob- 
served Hi column density would thus be Ahi ~ 
several 10^'^ cm^^. High sensitivity Hi 21cm 
observations of Virgo spiral galaxies with a dis- 
torted atomic gas distribution (see e.g. Phookun 
& Mundy 1995) have shown that low surface 
brightness structures appear when the sensitivity 
reaches several 10^^ cm~^. This could be also the 
case for mildly Hi deficient galaxies DEF < 0.4, 
where re- accretion is important (see Table 1). 



We suggest that most of the stripped gas, which 
can be found beyond the galactic disk, is thus very 
hkely either hot (T > 10® K), ionized, and has a 
fow density (n > rijcM), or molecular, cold (T ^ 
10— 100 K) and has a high density {n > 100 cm^"^). 
In the second case, only ~10% of the total gas 
column density is in atomic form due to X-ray 
dissociation. 

Combined Hi, CO, Ha, and X-ray observations 
will give informations about the different phases 
of the ISM. High resolution interferometric ob- 
servations of single cluster galaxies are necessary 
because they can be directly compared with our 
model snapshots. The detailed comparison be- 
tween simulations and observations will be a cru- 
cial test for our model. It will allow to draw con- 
clusions about the possible phase transitions (see 
above) of the ISM during the ICM-ISM interac- 
tion. We want to stress here that it is very im- 
portant to compare both, the gas distribution and 
the velocity field in order to have as much con- 
straints as possible. First results were already ob- 
tained with this technique: For the Virgo clus- 
ter sphal galaxy NGC 4522 VoUmer et al. (2000) 
succeeded in explaining extra planar Ha emission 
with a ram pressure scenario using one of the sim- 
ulation snapshots. VoUmer et al. (2001) compared 
multi wavelength observations of the Coma clus- 
ter galaxy NGC 4848 with a simulation snapshot. 
They showed that all observations are consistent 
with the ram pressure scenario. Studies of other 
cluster galaxies are in progress. The observations 
of single cluster galaxies will show the way how to 
model the properties of the ISM (e.g. phase transi- 
tions, ionization, temperature) and their evolution 
during an ISM-ICM interaction. 

6. Results 

For the analysis of the simulation we concen- 
trate on two time dependent parameters which 
are of particular relevance for the evolution of the 
galaxy: 

• the total gas mass Mtot within the galaxy's 
radius i?=20 kpc and a constant height H^l 
kpc, 

• and the central Hi surface density of the 
galaxy Sccntrai with a spatial resolution of 
^2 kpc. This is a slightly different definition 



than that of Cayatte et al. (1994) (Table 1), 
which loses its meaning if D^n < 1525 /4. 

The time dependence of these quantities is plotted 
in Fig. 13. For clarity, only the curves correspond- 
ing to Pram ~ 10 and lOOpo'^^o ^^'^ shown. The 
other curves are located between them. The differ- 
ent inclinations angle Q between the disk and the 
orbital plane are plotted with different lincstyles. 
The maximum ram pressure strength of lOOpo^o 
corresponds to the lower thick graphs, the strength 
of lOpo^o corresponds to the upper graphs in the 
same panel. Wc recall here that in these simula- 
tions the maximum ram pressure (i.e. the galaxy's 
closest approach to the cluster center) occurs at 
t — yr, indicated by the vertical line in Fig. 13. 

6.1. The total gas mass 

Observationally this is the most important pa- 
rameter because it gives directly the Hi deficiency, 
which is defined as 



DEF =< log X >T,D - log X. 



obs 



(18) 



(GiovancUi & Haynes 1985), where Xobs is the 
observed quantity of the cluster galaxy and < 
log X >T,D is averaged over a sample of field 
galaxies of morphological type T and optical di- 
ameter D. X can be the Hi mass, the Hi mass 
divided by the luminosity, or the Hi mass divided 
by the square of the optical linear diameter. In 
each case X{T, D) is proportional to the Hi mass 
of a galaxy. Since in our model the galaxy's di- 
ameter and luminosity do not change during the 
simulations, we define the model Hi deficiency as 
follows: 



DEF 



log(M^V<i) 



(19) 



where M^j and M^j are the total Hi masses at the 
beginning and the end of the simulation, i.e. after 
t=500 Myr. We want to point out three important 
issues: 

(i) In the model, ram pressure stripping acts 
very rapidly. For a cloud at rest, the timescale to 
be accelerated to the escape velocity Wcsc is given 
by Murakami & Babul (1999): 



VcscMcX 



Pr. 



.^Rt 



(20) 



where Md is the cloud mass, i?ci is the cloud ra- 
dius, and Pram is the ram pressure of the hot ICM. 
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Assuming Wcsc = v2wKopior =200 kms ^, _Rci =10 
pc, Mci = lO'^ Mq, and Pram = lOOpo^o gives 
r ~ 3 X 10^ yr. This has also been pointed out 
by Abadi et al. (1999). However, for a Lorentzian 
temporal profile ram pressure acts over a long pe- 
riod (up to several 10* yr). Especially in the case 
of face-on stripping (8 ~ 90°), where the restoring 
forces due to the galaxy's gravitational potential 
are minimum, stripping begins t ^ 300 Myr be- 
fore the closest passage to the cluster center. For 
decreasing this offset decreases. 

(ii) The consequences of ram pressure stripping 
on the total gas Hi content can be observed ear- 
lier for higher maximum ram pressure values and 
higher inclination angles (Fig. 13 and Table 1). 
For a maximum ram pressure of lOOpoWo stripping 
begins ^ 300 Myr before ram pressure reaches its 
maximum, whereas for a maximum ram pressure 
of 20/9oi'o and an inclination angle of 8 = 0° strip- 
ping begins ~50 Myr after ram pressure reaches its 
maximum. The effects of ram pressure stripping 
on non Hi deficient galaxies can thus be observed 
in some cases before the galaxy's closest passage 
to the cluster center. In the case of nearly edge- 
on stripping O < 30° the gas in the outer disk 
is mainly pushed to the inner part of the galactic 
disk in the beginning of the stripping (this cor- 
responds to the wings of the Lorentzian profile). 
In this stage of the ISM-ICM interaction the gas 
distribution in the galactic disk can be very asym- 
metric. 

The earliest moment when an Hi deficient 
galaxy {DEF > 0.3) with a distorted gas distribu- 
tion can be observed is t ^ —170 Myr (for face-on 
stripping) . If one assumes a constant mean veloc- 
ity of Vgai = 1700 kms"^ with respect to the ICM 
this corresponds to a minimum projected distance 
of ~ 1° from the cluster center. For 8 < 45° where 
there is a pronounced asymmetric gas distribution 
within the galactic disk t ^ —40 Myr which trans- 
lated to a minimum projected distance of ^ 15'. 

Fig. 14 resumes the situation. The critical dis- 
tance at which the total enclosed mass within 
20 kpc drops below 5 x 10^ Mq (solid lines) and 
3 X 10^ Mq (dashed lines) is plotted as a func- 
tion of the normalized ram pressure maximum for 
each simulation. Negative values indicate that the 
galaxy falls into the cluster center. The different 
symbols correspond to the different inclination an- 
gles 8 (triangles: 8 = 90°, diamonds: 8 = 45°, 



stars: 8 = 20°, crosses: 8 = 0°). The miss- 
ing points or lines in the upper part of the graph 
are due to the fact that the enclosed mass never 
drops below the given limit, i.e. ram pressure is 
not strong enough. The maximum projected dis- 
tance at which a distortion of the gas distribution 
of a non Hi deficient galaxy can be observed when 
it falls into the cluster core is ~ 2° (only in the case 
of face-on stripping). The maximum distance for 
inclination angles 8 < 45° or Hi deficient galaxies 
is - 1°. 

(iii) For low ram pressure maxima and low in- 
clination angles a considerable amount of gas ac- 
cretes back during 200-300 Myr (Fig. 13, Table 1). 
This accretion sets in at t ^100-200 Myr after the 
closest approach to the cluster center. In the case 
of edge-on stripping the accreted gas mass can at- 
tain values up to 5 X 10* M©. This corresponds to 
a maximum mass accretion rate of Af ~1 M© yr~^. 
The fraction of the stripped gas which falls back 
to the galaxy depends primarily on the inclination 
angle 8 between the disk and the orbital plane. 
For one given 8 this fraction increases with de- 
creasing maximum ram pressure. In the case of 
edge-on stripping more than 50% of the stripped 
gas can fall back to the galaxy, whereas in the 
case of face-on stripping only a few percent of the 
stripped gas re-accretes. The absolute value of the 
gas which falls back onto the galaxy as a function 
of the maximum ram pressure and the inclination 
angel 8 can be approximated crudely: 



Maccr ^ 6X ((-^)„,axXsin(4(0+lO°))') 

V pcivi 10 / 



-0.7 



(21) 

Thus, the re-accreted mass is approximately pro- 
portional to (pz;^)~°'^, where v± is the component 
of the galaxy's velocity within the cluster perpen- 
dicular to the disk plane. 

6.2. The central surface density 

If one assumes a Schmidt law (Schmidt 1959, 
Kcnnicutt 1983) for the global star formation rate 
of a galaxy, the central surface density Sc is di- 
rectly linked to the star formation activity in the 
central part of the galactic disk. Due to the dis- 
creteness of the model and the missing resolution 
in the galaxy center this parameter has a large 
scatter. For all simulations Sc stays constant at 
the initial value until t ^ —200 Myr. In the case of 
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face-on stripping (6 = 90°) Ec begins to decrease 
at i ~ -100 Myr. It drops within At ~ 150 Myr to 
its minimum value and then stays constant. The 
final central surface density decreases with increas- 
ing maximum ram pressure. For the highest ram 
pressure maximum (pram = lOOpo^o) i^ drops to 
about one third of its initial value. 

For e = 45° Sc begins to rise at t ~ -100 Myr, 
because a part of the gas is pushed to smaller 
galactic radii in the beginning of the stripping. 
The maximum which is about 1.2 times the initial 

value is reached at t ~ —100 50 Myr. For the 

lowest ram pressure maximum (pram = lOpo^o) 
there is a second increase and a second maxi- 
mum (1.4 times the initial value) at i ~ 230 Myr. 
For higher ram pressure maxima Sc decreases for 
t > Myr. 

For 6 = 20° Ec begins to rise later but more 
rapidly than for 6 = 45° (at t ^ -100 Myr) when 
the gas is pushed to smaller galactic radii. The 
maximum is about 1.4 times the initial value. 
Only for the highest ram pressure maximum 
(Pram — lOOpo'^^o) ^c drops bclow its initial value 
for t > Myr. In the simulations with lower ram 
pressure maxima Ec has a secondary maximum at 
t ~ 250 Myr which has the same value than the 
first one. 

In the case of edge-on stripping (9 = 0°) Sc 
begins to rise rapidly at i ^ —50 Myr and reaches 
its maximum shortly after i = Myr. The time 
when it reaches the maximum increases with de- 
creasing ram pressure maximum. The maximum 
of Ec increases with increasing ram pressure max- 
imum Pram- For Pram = 100 po^'o i* rcachcs 1.6 
times its initial value. As for Q = 20° there is a 
secondary maximum at t ~ 320 Myr except for 
Pram = WO pqVq . The offsct between the two max- 
ima is constant At ^ 300 Myr, i.e. the secondary 
maximum appears later for decreasing Pram- 

In order to summarize our findings, we conclude 
that the central surface density Sc decreases af- 
ter the galaxy's passage near the cluster center in 
the case of high inclination angles (nearly face-on 
stripping). For small inclination angles O (nearly 
edge-on stripping) Ec rises and reaches a maxi- 
mum when the galaxy passes the cluster center. 
In addition, there can be a secondary maximum 
after At - 300 Myr. 



7. Discussion 

7.1. Comparison with Hi observations 

The results of our simulations can be directly 
compared with the Hi 21 cm observations of Virgo 
cluster spirals made by Cayatte et al. (1994). We 
calculated the Hi diameter in the same way as 
these authors: it is defined by a limiting Hi col- 
umn density of 10^° cni^^. This parameter is plot- 
ted in Fig. 15. It shows the normalized Hi to the 
optical diameter of our model together with the 
observed values as a function of the Hi deficiency. 
The model values are taken from the final states 
of our simulations. We assumed an optical diam- 
eter of 20 kpc (which corresponds to ^7.5 disk 
scale lengths) for our model galaxy. The initial Hi 
diameter is Dgj'=30 kpc. 

There is an excellent agreement between the 
model values and the observational data for the 
normalized Hi to optical diameter. We are able to 
reproduce the slope and the scatter of the observed 
data. 

Now, if we assume that the Hi surface den- 
sity has the same value before and after the strip- 
ping event, the final Hi mass is given by M^" ~ 
7rEHi(-RH")^, where i?^" is the final Hi radius at 
the end of a simulation. Thus the Hi deficiency 
depends on the final Hi diameter in the following 
way: DEF ex log(-Dj^j ). This curve is represented 
by the solid line in Fig. 15. If we assume that the 
Hi surface density is proportional to R~^ and does 
not change during the stripping event, we obtain 
DEF oc log(I?jjj^). This curve is represented by 
the dashed line in Fig. 15. Thus, the modeled and 
observed Hi diameters as functions of the Hi de- 
ficiency are consistent with a constant radial Hi 
surface density, which has, in a first approxima- 
tion, the same value before and after the stripping 
event. However, Fig. 13 shows that the final sur- 
face density can be decreased by up to a factor of 
^3 with respect to the initial surface density. 

Cayatte et al. (1994) divided the observed spi- 
ral galaxies into four groups: 

Group I: galaxies with Dm/Do > 1.3, where 
Dm and Dq are the Hi and optical diameters re- 
spectively. These are cases of rather unperturbed 
Hi disks. 

Group II: galaxies with 0.75< Dyh/Dq < 1.3. 
The surface density falls off faster beyond the op- 
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tical half radius. 

Group III: galaxies with Dhi/^'o < 0.75. The 
disks are strongly truncated. 

Group IV: anemic galaxies with a central Hi 
hole. 

Since the model is in principle not able to re- 
produce a central hole which might be due to 
mechanisms other than ram pressure stripping, 
we will not discuss group IV. Anemic galaxies will 
be discussed in Section 8.2. The model surface 
densities were fitted with a fourth order polyno- 
mial to obtain smooth curves. This procedure 
can give overestimated values of the central sur- 
face density which are therefore not comparable 
to the value determined in Section 6.2. However, 
the overall behavior of the profile (especially the 
cut off) is well reproduced. We divided the final 
model surface density profiles by the initial one. 
These graphs are shown in Fig. 16. These normal- 
ized profiles can be directly compared with those 
of Fig. 6 in Cayatte et al. (1994). We can now 
assign each model density profile to one of the 
three groups. This gives: 

Group I: RUN A(10,0), B(20,0), C(50,0), D(100,0), 

E(10,20), F(20,20) 

Group H: RUN G(50,20), H(100,20), 1(10,45), 

J(20,45), K(50,45), M(10,90) 

Group HI: RUN L(100,45), N(20,90), 0(50,90), 

P(100,90) 

The membership of a simulation to a group de- 
pends on the form of the orbit, i.e. on the max- 
imum ram pressure and on the inclination angle 
between the galaxy's plane and the orbital plane. 
The comparison with the observed data shows that 
we can reproduce the scatter of the data over the 
whole range of deficiencies. The absolute values 
of the surface density profiles for the group I and 
III galaxies can be reproduced, whereas the values 
for the group II (intermediate deficiency) galax- 
ies are higher than the observed values. Only 
one simulation (H(100,20)) of a group II galaxy 
shows a significantly decreased final surface den- 
sity (Shi < 10^^ cm~^). Thus, our data does not 
show a systematic decrease of the central Hi sur- 
face density after the stripping process for group 
II galaxy as observed by Cayatte et al. (1994). We 
suggest that this decrease is related to the physical 
environment in the inner stellar and gaseous disk 



as the UV radiation field, the supernova rate, or 
the metalicity of the gas which influence the ex- 
tent of the photodissociation regions and the gas 
heating and cooling rates leading to different Hi 
to H2 ratios. 

7.2. Comparison v^^ith SPH simulations 

Abadi et al. (1999) made ram pressure strip- 
ping simulations using an SPH code. They mod- 
eled a spiral galaxy moving with a constant ve- 
locity within an ICM of constant density. This 
corresponds to a circular galaxy orbit within the 
cluster. Therefore they do not have re-accretion 
of material which is mainly due to the decrease of 
ram pressure when the galaxy leaves the cluster 
core. Their timescale for ram pressure stripping 
is between 1 and 5x10^ yr, thus comparable to 
our simulations. In difference to our model they 
have used a decreasing gas surface density profile 
Sgas oc R~^ that corresponds to the total gas sur- 
face density (atomic and molecular). They there- 
fore underestimated the final Hi deficiency. This 
becomes clear if one estimates the Hi deficiency us- 
ing their Hi to optical radius fraction and Fig. 15. 
For example, i?Hi/-Ropt ~ 0.3 for their run D and 
E. According to Fig. 15 the resulting Hi deficiency 
is DEF > 0.8 compared to the given value of 
DEF - 0.4. 

7.3. Interpretation of the Hi data 

We conclude that our simulations can repro- 
duce the observational data of Cayatte et al. 
(1994) in its behavior, absolute values, and scat- 
ter. Therefore, we suggest that ram pressure strip- 
ping is the main cause of their Hi deficiency. 

It is worth noting that we can reproduce the Hi 
observations using only one passage of the galaxy 
near the cluster center. As a typical orbital pe- 
riod of an Hi deficient spiral galaxy is of the order 
of several Gyr, it is possible that it has already 
passed through the cluster center more than once. 
The galaxy trajectory model has shown that the 
existence of a perturbing gravitational potential 
(M86) in the cluster center has two effects: (i) It 
can make the central mass concentration (M87) 
move. As the intracluster gas will follow the vari- 
able gravitational potential, the location of high- 
est gas density within the cluster can also change, 
(ii) The galaxy orbits will change leading to a ram 
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pressure profile changing with time for each pas- 
sage in the cluster center. 

Our simulations indicate (Fig. 15) that the 
lower Hi surface densities of galaxies with DEF ^ 
0.5 is due to internal parameters of the galaxy 
(Section 7.1) which lower the Hi to H2 ratio within 
the optical disk, whereas the total gas surface den- 
sity stays the same. Furthermore, the simulations 
of Valluri (1993) indicate that the gravitational 
potential of the galaxy does not change signifi- 
cantly when the galaxy passes near the cluster 
center. For the estimation of the final Hi defi- 
ciency after a stripping event, the stellar and gas 
surface densities can thus be regarded as constant 
in a first approximation. A galaxy with an al- 
ready truncated Hi radius with an unchanged to- 
tal gas surface density will only lose more gas in a 
further stripping event if the maximum ram pres- 
sure and/or its inclination angle with respect to 
the orbital plane is higher than before. If a galaxy 
has already accomplished several orbits within the 
cluster, the orbit leading to the closest approach to 
the cluster center and/or having the highest incli- 
nation angle O is decisive for its Hi deficiency. On 
the other hand, if the final surface density is de- 
creased, a further passage through the cluster core 
will result in a higher Hi deficiency. The observed 
population of deficient galaxies represents galax- 
ies which have passed through the cluster center 
at least once. 

The different groups of Cayatte et al. (1994) 
represent galaxies on different orbits in the cluster. 
The fact that we can assign galaxies in the final 
state after the interaction to these groups confirms 
the conclusion that all Hi deficient galaxies have 
already passed through the cluster center. 

The effects of ram pressure can be observed be- 
fore the closest passage of the galaxy to the clus- 
ter center only in the case of inclination angles 
G > 45°. The observed asymmetries of the gas dis- 
tribution of Hi deficient {DEF > 0.3) in the Virgo 
cluster which are located more than 1° from the 
cluster center are due to a former ISM-ICM inter- 
action. This interaction has already ceased but its 
consequences are still observable. This is also the 
case for non Hi deficient galaxies which are located 
more than 2° from the cluster center. Thus, these 
galaxies showing important distortions in their Hi 
distribution have already passed through the clus- 
ter center. They are not infalling galaxies. 



7.4. The amount of stripping 

The knowledge of the dependence of the Hi de- 
ficiency on the maximum ram pressure and the 
inclination angle together with the galaxy's ra- 
dial velocity can constrain the galaxy orbits nec- 
essary to obtain a given Hi deficiency. We there- 
fore attempt to give a simple analytical formula. 
Fig. 17 shows the fraction between the final and 
initial total Hi mass as a function of the quantities 
(p v'^) / {po Vq) (ram pressure strength) and (in- 
clination angle). The solid line represents a linear 
fit. The least square fit for the fraction between 
the initial and final total Hi mass is 

0.25 (^^) sin2(A(e + io°)) + 0.84 

(22) 
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7.5. Galaxies with different properties 

The present simulations were made using a 
galaxy model with a constant rotation curve of 
''^rot = 140 kms^^, a total atomic gas mass of 
M^f = 6 X 10^ Mq, and an initial Hi diame- 
ter of Dhi=30 kpc. For a galaxy with different 
Wrot, -^m*' ^^"i ^Hi the resulting Hi deficiency 
changes in the following way: 

The final radius of the Hi gas distribution at 
the end of a simulation is given by 



-RhI — ^gasPram^: 
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(23) 



where Pram = PiCMWgaiaxy ^nd Sgas is the total gas 
surface density. The fraction between the initial 
and final Hi mass is given by 
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where _Ryj " is the initial/final Hi radius. The 
total gas surface density can be regarded as con- 
stant (Section 7.1) Sgas = S™*- With Mm = 
(i?|5?i*)2 one obtains 
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For a maximum ram pressure Pram, the fraction 
between the initial and final Hi mass for a galaxy 
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is the mass fraction for Ur, 



=140 km s" 



where (^^ 

M*if = 6 X 10^ Mq, and Dhi==30 kpc. In order to 
recover the deficiencies, one simply has to apply 
its definition Eq. 19. We made two simulations 
with Wrot=250 kms-\ M^f = 8 x 10^ M©, and 
£)hi=54 kpc: (i) fw^ = 20/3OWO and 6 = 20°, (ii) 



pv 



20po'i'o and O = 45°. In order to verify 



Eq. 26, we calculated ( 



M" 



Mil 



using Eq. 26 for 



these two simulations. The two resulting values 
can be seen as open triangles in Fig. 17. They are 
in good agreement with the mass fractions of the 
simulations of Tabic 1 . 

7.6. Possible local star formation mecha- 
nisms 

The most important consequence of ram pres- 
sure stripping on cluster spiral galaxies is the re- 
moval of a large fraction of their gas reservoir. 
Without the supply of fresh gas, the star formation 
activity of a stripped galaxy will decrease on the 
timescale of a cluster crossing time. Nevertheless, 
the interaction of the galaxy's gas with the ICM 
can lead to a short enhancement of its star forma- 
tion activity. In this section we investigate possi- 
ble star formation mechanisms triggered by a ram 
pressure stripping event. We first resume the rel- 
evant results of our simulations before discussing 
their consequences on star formation activity. 

• A considerable fraction of the initial galaxy's 
gas mass can be accreted back on the galaxy 
in the case of low ram pressure maxima and 
low inclination angles (up to 10%). 

• There is a strong enhancement of the central 
surface density for a high and intermediate 
maximum ram pressure and small inclina- 
tion angles O. 

• In case of edge-on stripping and high maxi- 
mum ram pressure, the central surface den- 
sity increases by a factor 1.5 during the 
galaxy's passage in the cluster center. 

If we relate the enhancement of the gas surface 
density to a local star formation rate via the 
Schmidt law (Schmidt 1959, Kennicutt 1983) and 
via star formation due to colliding flows (sec e.g. 
Hunter et al. 1986), these results lead to two pos- 
sibilities of an enhancement of star formation: 



(i) An instantaneous burst in the inner galactic 
disk which happens before {At < 10^ yr) or while 
the galaxy passes through the cluster center. For 
the Virgo cluster this happens within a radius of 
^ 40' around the cluster center. The enhance- 
ment of the star formation activity is strongest if 
the galaxy passes very close to the cluster center 
(~100 kpc) and if the stripping is edge-on. Using a 
Schmidt-law, we estimate that the star formation 
rate within a radius R = 0.25 x i?25 can increase 
up to a factor ^2. 

(ii) The re-accretion of gas leads to star forma- 
tion when the infalling gas clouds collide with the 
clouds in the disk. This can lead to a star forma- 
tion event within the whole galactic disk which be- 
gins ^ 2 X 10* yr and ends up to ^ 5 x 10* yr after 
the galaxy's closest passage to the cluster center. 
In the case of maximum re- accretion (small 0, low 
maximum ram pressure) ~ 10% of the total gas 
content falls back to the galaxy within ~ 3x 10* yr. 
If ^10% of this gas is turned into stars, the result- 
ing star formation rate is SFR ^0.2 M0 yr^^. 

We conclude that ram pressure stripping can in- 
duce two types of star formation activity: a short 
burst during its passage in the cluster center or a 
longer lived more continuous star formation activ- 
ity due to the re-accretion of the stripped gas. The 
strength of the star formation activity is governed 
by the galactic orbits in the cluster. This confirms 
the suggestion of Bothun & Dressier (1986) that 
ram pressure causes the Hi deficiency and can in- 
duce star formation. Several individual cases will 
be discussed in forthcoming papers. 

8. The influence of orbits on galaxy prop- 
erties 

With the acquired knowledge of the previous 
section we can now study the link between the 
properties of Virgo cluster spiral galaxies and their 
orbits. Our model leads to the following scenario: 
galaxies on radial orbits will be Hi deficient re- 
gardless of their morphology. These galaxies have 
already passed through the cluster center at least 
once. For most of the group I and II galaxies, 
which have entered the cluster a few crossing times 
ago, the passage with the smallest impact param- 
eter was decisive for their Hi deficiency. 

Dressier (1986) and more recently Solanes et al. 
(2000) showed that the velocity dispersion of Hi 
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deficient cluster galaxies decreases with increasing 
distance to the cluster center, indicating that these 
galaxies are on highly eccentric orbits. In addition, 
they found indications for a segregation among the 
orbits of infalling spirals according to their Hi de- 
ficiency regardless of galaxy morphology. Interest- 
ingly, Solanes et al. (2000) pointed out that the 
velocity dispersion of Hi deficient galaxies in the 
Virgo cluster does not decrease with radius. They 
interpreted this result as indicative for a gravita- 
tional field, which is strongly perturbed by ongo- 
ing mergers of major subclumps. 

In order to investigate the link between Hi defi- 
ciency and galaxy orbits, we use position- velocity 
plots of different galaxy populations in the Virgo 
cluster together with our model orbits. 

8.1. Hi deficient galaxies 

We show the Hi deficiency data collected by 
Guiderdoni & Rocca-Volmerange (1986) in a ra- 
dial velocity-projected distance plot (Fig. 18). 
The curves of the model orbits (Section 2) are 
plotted as dotted lines, the Keplerian and escape 
velocity for the gravitational potential of M87 as 
solid and dashed lines. The model orbits cover 
about half of the total area in the radial velocity- 
projected distance plane. The size of the crosses 
is proportional to the Hi deficiency. We only plot 
galaxies with DEF > 0.48. 9 out of 41 Hi defi- 
cient galaxies (~ 20%) lie beyond the model or- 
bits marked by the dotted lines for a maximum 
projected distance of 12° (NGC 4193, NGC 4212, 
NGC 4216, NGC 4235, NGC 4260, NGC 4419, 
NGC 4450, NGC 4522, NGC 4826). The range of 
the projected distance of these galaxies is between 
2° and 7°. Solanes et al. (2000) did not find a 
decrease of the velocity dispersion for distances 
smaller than 5° because of these high velocity Hi 
deficient galaxies. Their velocities are even higher 
than the Keplerian velocity due to the M87 grav- 
itational potential but still below the escape ve- 
locity. Thus it is likely that they evolve on bound 
orbits. We will discuss this issue in Section 9. For 
larger radii we find a trend that Hi deficient galax- 
ies evolve on radial orbits as already pointed out 
by Dressier (1986) and Giraud (1986). 



8.2. Anemic galaxies 

Anemic galaxies are a special class of galaxies 
with a very low arm inter-arm contrast (van den 
Bergh 1976). They are most similar to the passive 
spirals observed in high redshift clusters (see e.g. 
Poggianti et al. 1999) and are believed to repre- 
sent the oldest spiral galaxy population which has 
entered the cluster several Gyr ago. We suggest 
the following evolutionary scenario for these galax- 
ies: after a major stripping event the galaxy has 
lost its gas in the outer disk. The ongoing star for- 
mation activity exhausted its remaining gas reser- 
voir in its entire disk. Without the possibility of a 
supply of fresh gas from the outer disk, star forma- 
tion activity in the spiral arm decreased which led 
to a decreasing arm inter-arm contrast, i.e. the 
galaxy became anemic. Thus, an anemic galaxy 
has entered the cluster more than one crossing 
time ago on a radial orbit. In order to investi- 
gate if anemic galaxies are still on radial orbits or 
if tidal interactions have lead to more circular or- 
bits, we show the distribution of this population 
in the radial velocity-projected distance plot in 
Fig. 19. We include here also the group IV galax- 
ies of Cayatte et al. (1994), because we assume 
that the reduced Hi diameter regardless of a cen- 
tral hole is due to the major stripping event. 

The small number of orbits prevents us from a 
quantitative tracing of the dispersion velocity as a 
function of the projected distance. Nevertheless, 
we observe that 12 out of 15 anemic galaxies lie on 
radial orbits except (NGC 4192, NGC 4450, and 
NGC 4419). NGC 4192 is seen nearly edge-on. 
Therefore, it is difficult to estimate its arm inter- 
arm contrast. The low Hi deficiency of NGC 4192 
indicates that this galaxy might be misclassified 
as an anemic galaxy. If we exclude this galaxy, 
there is an indication that the velocity dispersion 
decreases with increasing distance. This is consis- 
tent with anemic galaxies being preferentially on 
radial orbits. As they represent a spiral popula- 
tion which entered the cluster more than several 
Gyr ago, there is an indication that their orbits 
were not perturbed by the changing gravitational 
field due to the merging of subclumps. If this ten- 
dency turns out to be real, this would mean that 
only the infalling galaxies are scattered with high 
velocities to large distances from the cluster cen- 
ter during the merging of a major subclump (e.g. 
M86) and not the old spiral population. Interest- 
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ingly, Solanes et al. (2000) showed that the veloc- 
ity dispersion of early type spiral galaxies does not 
decrease. Anemic galaxies might thus form a spe- 
cial class of galaxies whose orbits are of different 
shape than those of early type spirals. 

9. Infall 

We will now try to explain the existence of 
the group of high velocity Hi deficient galaxies 
which can be found at projected distances of ^ 4° 
(NGC 4193, NGC 4212, NGC 4216, NGC 4235, 
NGC 4260, NGC 4419, NGC 4450, NGC 4522). 
One of its members (NGC 4522) has been observed 
in Ha by Kenncy & Koopmann (1999). They 
claim that there are extraplanar ionized filaments 
and identify it as one of the best candidates for 
ram pressure stripping. As we have shown be- 
fore, spiral galaxies on radial orbits which lead to 
a strong stripping event do not permit these po- 
sitions in the radial velocity - projected distance 
plot. Thus, there must be another explanation. 

The basic idea comes from the presence of M86 
in the vicinity of M87. Since it has a high radial 
velocity with respect to M87, it is probable that 
M86 is falling into the cluster from behind. In 
addition, the central dwarf population shows al- 
ready a spatial distortion in the direction of M86 
(Schindler et al. 1999). We interpret this as the 
consequence of the passage of M86 near the cluster 
center, i.e. M87. If the subclump of M86 has en- 
tered the Virgo cluster with its own spiral galaxy 
population, these could have been captured by the 
Virgo cluster. The outcoming orbits of this in- 
falling spiral population can be very different from 
those which are already settled. 

We have tested this hypothesis with the help 
of our orbit simulations. M87 has been placed at 
rest in the center. M86 is located at (0.5 Mpc, 
Mpc, 5 Mpc) and with an initial velocity vec- 
tor (0 kms~^, 1 kms""'^, kms""'^). For the test 
galaxy we adopted: x = (0.5 Mpc, 1 Mpc, 5 Mpc) 
and V — (258 kms~^, 1 kms^^, 258 kms~^), i.e. 
the test galaxy is orbiting at a distance of 1 Mpc 
around M86 which is falling radially (along the 
LOS) into the gravitational potential of M87. We 
have included the full symmetric gravitational po- 
tential of both galaxies. The effects of an internal 
perturbation of M86 due to the tidal interaction 
with M87 is neglected. The assumption is justi- 



fied, because the difference between their radial 
velocities Aw ~ 1500 kms^^ is much larger than 
the dispersion velocity of M86. 

The resulting trajectories can be seen in Fig. 20. 
When the M86 clump approaches M87 the latter 
also begins to move. The test galaxy is orbiting in 
the joint gravitational potential of M87 and M86 
following M86 in the z-direction. When passing 
near M87 it can even be ejected from the cluster 
as it happens at the end of our simulation. The 
corresponding radial velocity - projected distance 
plot can be seen in Fig. 21. The positions and 
velocities are plotted with respect to the positions 
and velocities of M87, the center of the Virgo clus- 
ter. Clearly there is the possibility to obtain high 
radial velocities at large projected distances espe- 
cially when the galaxy is ejected from the cluster. 
There is also the effect that M87 begins to move, 
changing the position of the cluster center. 

We suggest that the existence of the special 
galaxy population with high radial velocities at 
large projected distances is linked to the trajectory 
of M86 perturbing the gravitational potential. In 
this scenario these galaxies pass the cluster cen- 
ter with high velocities. The increased velocity 
leads to an enhanced ram pressure (pram c>c vt^i)- 
Solanes et al. (2000) also suggested that gravita- 
tional perturbations due to major subclumps (as 
M86) are responsible for the high velocities Hi de- 
ficient galaxies. If these galaxies are bound, they 
will evolve on extremely eccentric orbits with very 
long periods. After the merging of the subclump, 
they will be most probably found at very large 
distances from the cluster center and the velocity 
dispersion of the remaining Hi deficient galaxies 
will again decrease. 

Since there are indications that galaxies are ac- 
creted in subclumps, the accretion rate is expected 
to be highly irregular showing strong peaks when a 
subclump as M86 is falling into the cluster adding 
new spiral galaxies to the existing population. 

10. Galaxy evolution within the cluster 

Accreted spiral galaxies are preferentially on ra- 
dial orbits which lead them near the cluster center 
where ram pressure is most efficient. Very recently 
accreted galaxies (~1 Gyr) can have an enhanced 
star formation rate in the outer disk due to the 
re-accretion of stripped gas. During the ram pres- 
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sure stripping event, the atomic gas is stripped 
from the outer disk leaving the molecular gas un- 
affected. Consequently the star formation activity 
in the outer disk region of highly stripped galaxies 
decreases, whereas the star formation rate in the 
inner disk stays roughly constant after a possible 
short increase (see Section 7.6). This leads to a 
decrease of the overall star formation rate. Koop- 
mann & Kenney (1998) noticed that the morpho- 
logical classification of the spiral galaxies in the 
Virgo cluster taking only into account the overall 
star formation rate is often ambiguous. Galaxies 
classified in this way as early type spirals become 
late types if one takes into account their disk to 
bulge ratios in the R band. These galaxies show 
a reduced overall star formation rate, which is 
mainly due to a strongly reduced star formation 
rate in the outer disk. We suggest that this spe- 
cial type of galaxies has undergone a ram pressure 
stripping event. 

We will now estimate the influence of ram 
pressure stripping on the global star formation 
rate. We assume a gas surface density distribution 
Sfja = 5]o exp(— i?/a) for the molecular phase and 
Efji = const for the atomic phase within an opti- 
cal radius. Furthermore, the star formation rate 
is given by a Schmidt law (Kennicutt 1998) 

SFR^ 2^2.5x10-^" Rll\;^^dR, (27) 

Jo 

where Stot = Shs + Shi- It is assumed that star 
formation takes only place within the optical ra- 
dius. After a major stripping event, the Hi is 
-Rhi = 0.3 X i?opt. We consider two different cases: 

(i) So = 100 Mopc-^ a ^2 kpc, I;hi=5 Mqpc-^. 
For i?i = i?25 the total molecular gas mass and 
the total atomic gas mass within the optical radius 
are comparable Mh2 — M^a- The star formation 
drops by a factor 2 due to ram pressure stripping 
(from 3.6 to 1.7 Mgyr^^). 

(ii) So = 50 M0PC-2, a = 2 kpc, Shi=5 IVIgpc'^. 
In this case M-^^ ~ 0.5 x Mhi. The star formation 
drops almost by a factor 3 due to ram pressure 
stripping (from 2.2 to 0.8 Mgyr"^). 

For this high maximum ram pressure star for- 
mation induced by re-accretion is negligible. In 
the case of a central starburst triggered by ram 
pressure, we assume an increase of the local star 
formation rate within R < 0.25 i?25 by a factor 4 
during 5 x 10^ yr (Fig. 13. This leads to a mild 



increase of the global star formation rate (up to 
a factor 2)). The consumed molecular gas masses 
for both cases are (i) Mcons ^ 3 x 10^ Mq and 
(ii) Mcons ^ 1-3 X 10^ Mq. This corresponds to 
roughly 10% of the total molecular gas content in 
both cases. Thus, the occurrence of a central star 
burst can additionally lower the global star forma- 
tion rate by ~10% after a strong stripping event. 

If the gas content within an optical radius of a 
non deficient galaxy is dominated by the atomic 
phase, ram pressure stripping can account for the 
quenching of star formation at a timescale of ^1 
Gyr described by Poggianti et al. (1999). In the 
case of a strong edge-on stripping event, the global 
star formation will increase by a factor ~1.5 within 
^ 10^ yr and drop then by a factor 2-3 three 
with respect to its initial value. On the contrary, 
ram pressure cannot account for the morpholog- 
ical transformation, which has a larger timescale 
(several Gyr). This transformation might be due 
to the combined effect of high speed encounters 
with massive galaxies and the interaction with the 
gravitational potential of the cluster ( "galaxy ha- 
rassment"; Moore et al. 1996, 1998). 

We will now attempt to determine the galaxies' 
location in the z direction with respect to the clus- 
ter center (M87) using the following results from 
this article: 

• If an Hi deficient galaxy shows enhanced star 
formation in the disk due to re-accretion of 
the stripped material, it has been accreted 
recently (less than 1 Gyr) by the cluster. 
We observe the galaxy when it has already 
passed the cluster center. It is in the phase of 
re-accretion of the stripped material which 
collides with the galaxy's ISM leading to an 
enhancement of its star formation activity. 

• If an Hi deficient galaxy shows a perturbed 
Hi velocity field or a distorted Hi distribu- 
tion, it has passed the cluster center recently 
and is now coining out of the cluster center. 

• If a non Hi deficient galaxy shows a per- 
turbed Hi velocity field or a distorted Hi dis- 
tribution and is located more than 2° away 
from the cluster center, it has passed the 
cluster center recently and is now coming out 
of the cluster center. 
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These galaxies have in common that they are com- 
ing from the cluster core with the closest approach 
to the cluster center less than 1 Gyr ago. The sign 
of their radial velocity gives us the clue to place 
them with respect to M87. If the galaxy's radial 
velocity is positive it is located behind the cluster 
center, if it is negative it is located in front. 

11. Summary and conclusions 

In this work we have made N-body simulations 
in order to simulate the neutral gas content of spi- 
ral galaxies entering a cluster. We are interested 
in the effect of the ram pressure exerted by the 
hot intraclustcr medium on the ISM of the fast 
moving galaxy. At the same time we simulated 
the gas dynamics of a spiral galaxy on a radial or- 
bits within the gravitational potential of the Virgo 
cluster. We investigated different orbits varying 
systematically the inclination angle 8 between the 
disk and the orbital plane. The main results are: 

1. The amount of the stripping and thus the Hi 
deficiency depends on the galaxy orbit (min- 
imum distance to the cluster center, max- 
imum velocity, and inclination of the disk 
with respect to the orbital plane). We gener- 
alize the formula of Gunn & Gott (1972) in- 
cluding explicitly the inclination angle. We 
give an approximation of the relation be- 
tween the Hi deficiency on the one hand and 
the maximum ram pressure and inclination 
angle on the other hand. 

2. We have used a realistic temporal ram pres- 
sure profile, i.e. ram pressure decreases af- 
ter the closest passage to the cluster cen- 
ter. For maximum ram pressure of Pram — 
1000 - 5000 cm-3 (kms-i)^ and low incli- 
nation angles 8 < 45° we observe a phase 
of re-accretion where the gas, which has not 
been accelerated to the escape velocity, falls 
back to the galaxy. The fraction of the re- 
accreted mass can attain up to 10% of the 
total atomic gas mass. 

3. For maximum ram pressure of Pram > 
2000 cm^'^ (kms^^)^ and small inclination 
angles 8 < 20° the central gas surface den- 
sity can increase up to a factor 1.5 within a 
few lO^'yr. Assuming a Schmidt law for the 
star formation rate leads to an increase of 



the local star formation rate up to a factor 

-2. 

4. Deformations of the neutral gas content are 
mainly observable after the closest approach 
to the Virgo cluster center. Thus, the major- 
ity of galaxies with a peculiar Hi distribution 
are coming out of the cluster core. 

5. The observationally established correlation 
between the Hi to optical diameter and the 
Hi deficiency is well reproduced by our sim- 
ulations. 

We conclude that the scenario where the Hi defi- 
ciency is due to ram pressure stripping is consis- 
tent with all available Hi 21 cm observations of 
the Virgo spiral galaxy population. 
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Fig. 1. — The galaxy velocity as a function of 
the ICM density for a given ram pressure Pram = 
PiCMi'gaiaxy Solid line: the minimum ram pres- 
sure to affect the atomic gas disk Pram = 0.5 x 
10^ cm~^(kms~"'^)^. The area below the solid 
line corresponds to the part of the parameter 
space where ram pressure is not effective. Dot- 
ted line: Pram — 10^ cm~^(kms~^)^. Dashed line: 
Pram = 2 X fO"^ cm^"^(km s^^)^ . Dot-dashed line: 
Pram = 5 X 10^ cm^'^(km s^^)^ . Dot-dot-dashed 
line: Pram — 10'' cm^''(kms^^)^. 




Fig. 3. — Normalized ram pressure (in units of 
PoVq) as a function of time for the galaxy orbits 
shown in Fig. 2. 




" 


A 


2 


a . 





^^ 


k^^ 


2 


/ / 

'^ \ 


- 


^ 


: 



-2-10 1 2 
X-axis (Mpc) 



-2-10 1 2 
x-Qxis (Mpc) 



Fig. 2. — Test particle trajectory calculated with 
the fixed gravitational potential of M87 and M86. 
Left: a;-y-plane. Right: x-z-plane. This simu- 
lation does not trace the trajectory of one single 
galaxy but shows a variety of orbits corresponding 
to different initial conditions. 
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Fig. 4. — Ratio of the ram pressure and thermal 
pressure for the galaxy orbits shown in Fig. 2 as 
a function of the projected distance to the cluster 
center. 
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Fig. 5. — Linc-of-sight (LOS) velocity with re- 
spect to the cluster mean velocity as a function 
of the projected distance to the cluster center for 
the galaxy orbits shown in Fig. 2. 





distance (kpc) 

Fig. 6. — Decomposition of the rotation curve af- 
ter a disk evolution of 10^ yr. Dashed: halo, solid: 
bulge, dotted: disk; thick: total, solid with error 
bars: initial model rotation curve. The error bars 
represent the velocity dispersion of the N-body 
model. 



Fig. 8. — Illustration of the effect of ram pres- 
sure. The direction of the wind created by the 
galaxy's motion in the ICM is indicated by the ar- 
rows. Only clouds which are directly exposed to 
the wind are pushed by ram pressure (clouds no.l, 
2, 5, and 6). The dashed line represents the galac- 
tic plane, O is the inclination between the galactic 
plane and the wind direction. 
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Fig. 7. — Hi cloud distribution in z direction after 
a disk evolution 10^ yr. 
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Fig. 9. — Dashed: ram pressure profiles of Fig. 3. 
Solid: analytical ram pressure profile (Lorentzian, 
see Eq. 11). 
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Fig. 10. — Final Hi radius Rm = -Dhi/2 of our 
simulations as a function of the estimated radius 
i?GG using the formula of Gunn & Gott (1972). 
The solid fine corresponds to Rqg — Rm- 
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Fig. 11.— Snapshots of RUN C(50, 0). The 
elapsed time is indicated at the top of each panel. 
The galaxy is seen face-on and is moving in the 
south-east direction, i.e. the wind is coming from 
the south-east. This wind direction is indicated 
by the arrows whose length is proportional to pv'^ . 
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Fig. 12.— Snapshots of RUN J(20, 45). The 
elapsed time is indicated at the top of each panel. 
The galaxy is seen edge-on and is moving in the 
south-east direction, i.e. the wind is coming from 
the south-east. This wind direction is indicated 
by the arrows whose length is proportional to pv^ . 
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Fig. 13. — Variation of different galaxy parameters 
versus time. Upper panel: Total gas mass within 
20 kpc Mtot in M0. Lower panel: Central surface 
density Sccntrai in cm^^. In each panel two graphs 
are shown for the same inclination angles between 
the disk and the orbital plane corresponding to 
Pram = lOOpo^'o (thick huc) and p^am = lOpoi'o 
(thin line). The different inclination angles be- 
tween the disk and the orbital plane are: solid: 
90°; dotted: 45°; dashed: 20°; dash-dotted: 0°. 
The vertical line at t—0 yr indicates the closest 
passage of the galaxy to the cluster center. 
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Fig. 14. — Critical distance at which the total en- 
closed mass within 20 kpc drops below 5 x 10^ Mq 
(solid lines) and 3 x 10^ M0 (dashed lines) as a 
function of the normalized ram pressure maximum 
for each simulation. Negative values indicate that 
the galaxy falls into the cluster center. The differ- 
ent symbols correspond to the different inclination 
angles & (triangles: G = 90°, diamonds: 8 = 45°, 
stars: 9 = 20°, crosses: 9 = 0°). 





0.1 1 

0.0 



Fig. 15. — Normalized Hi to optical diameter as 
a function of the Hi deficiency. Squares: observed 
values (Cayatte et al. 1994); stars: model val- 
ues. The horizontal line represents the error in 
the Hi deficiency determination. The solid line 
corresponds to DEF ex log(£'jjj^), the dashed line 
corresponds to DEF oc log(I?jjj^). 
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Fig. 16. — Smoothed normalized model surface density profiles. The distance to the galaxy center is divided 
by the optical radius. The inclination angle between the galaxy's disk plane and its orbital plane together 
with the maximum ram pressure are indicated on top of each panel. The RUN is indicated in the lower right 
corner, the group membership (Cayatte et al. 1994) in the upper right corner. 



27 




pv'.sin'(9/lO.{0+lO°)) 

Fig. 17. — Hi deficiency as a function of tlie quan- 
tities (p v'^) I (pq vf^) (ram pressure strengtli) and 
Q (inclination angle). Tlic solid line represents a 
linear least square fit. The triangles correspond 
to simulations of galaxies with different initial Hi 
mass, initial Hi diameter, and rotation curve. 
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Fig. 19. — Radial velocity - projected distance plot 
for the anemic Virgo cluster spirals. The curves 
for for the model orbits (Section 2) are plotted as 
dotted lines. 
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Fig. 18. — Radial velocity - projected distance plot 
for the Hi deficient {DEF > 0.48) Virgo cluster 
spirals. The curves of the model orbits (Section 2) 
are plotted as dotted lines, the maximum Keple- 
rian velocity for the gravitational potential of M87 
as solid lines. The size of the crosses is propor- 
tional to the Hi deficiency. 




Fig. 20. — Trajectories of the test galaxy orbiting 
around M86 which is falling radially into the grav- 
itational potential of M87. Solid line: test galaxy. 
Dotted line: M87. Dashed fine: M86. 
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Fig. 21. — Radial velocity - projected distance plot 
for the trajectories of Fig. 20. The Keplerian ve- 
locity due to the gravitational potential of M87 is 
plotted as a thick line. The positions and veloci- 
ties are plotted with respect to the positions and 
velocities of M87. 
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Table 1 
Parameters and results of the different simulations 



RUN 


1' p«^ \a 
V po t^Q z max 


Qb 


At'' 


AD'^ 


Dh 


^central 


DBFS 


^"accr 


Maccr/M,tr ' 


A(10,0) 


10 





- 


- 


32 


12.71 


0.01 


2.4 


1.41 


B(20,0) 


20 





- 


- 


34 


13.85 


0.02 


3.1 


0.94 


C(50,0) 


50 





100 


180 


36 


15.31 


0.05 


5.7 


0.68 


D(100,0) 


100 





70 


125 


37 


12.76 


0.10 


5.1 


0.32 


E(10,20) 


10 


20 


50 


90 


33 


15.87 


0.08 


1.8 


0.17 


F(20,20) 


20 


20 








30 


16.26 


0.13 


2.5 


0.15 


G(50,20) 


50 


20 


-30 


-50 


20 


16.57 


0.23 


1.0 


0.04 


H(100,20) 


100 


20 


-200 


-360 


19 


4.9 


0.61 


0.4 


0.01 


1(10,45) 


10 


45 


-110 


-200 


20 


17.68 


0.34 


2.6 


0.08 


J(20,45) 


20 


45 


-130 


-230 


14 


18.03 


0.73 


1.5 


0.03 


K(50,45) 


50 


45 


-160 


-285 


14 


10.93 


0.72 


0.8 


0.02 


L(100,45) 


100 


45 


-170 


-300 


9 


4.83 


1.13 


0.3 


0.01 


M(10,90) 


10 


90 


-320 


-570 


14 


13.12 


0.60 


2.3 


0.05 


N(20,90) 


20 


90 


-320 


-570 


11 


10.75 


0.82 


1.2 


0.02 


0(50,90) 


50 


90 


-320 


-570 


8 


5.64 


1.11 


0.5 


0.01 


P(100,90) 


100 


90 


-320 


-570 


6 


3.33 


1.35 


0.3 


0.01 



^ Maximum ram pressure in units of po^o- 

^ Inclination angle with respect to the orbital plane in degrees. 

'^ Time between the beginning of the stripping (total gas mass within 20 kpc Mtot < 5 x 10^ M0) and the 
galaxy's closest passage to the cluster center in Myr. 

'^ Radial distance in kpc corresponding to At assuming a constant velocity of the galaxy with respect to 
the cluster center of Ugai=1700 kms^^. 

° Final Hi diameter in kpc. 

^ Final central surface density in Mq pc^^. 

s Final Hi deficiency. 

^ Re-accreted gas mass in 10^ ^'i-Q- 

' Re-accreted gas mass divided by the stripped gas mass. 
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